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1 Introduction 



The string fragmentation model |I[ is able to give a generally very good description of 
the distribution of hadrons in quark and gluon jets 0. The abundance of mesons with 
different flavour and spin can be described with a limited number of parameters, which have 
a natural interpretation. Thus the production of a qq pair can be understood as a tunneling 
phenomenon ||, which gives a suppression of strange quarks compared to the lighter u- and 
d-quarks. Also the relation between vector and pseudo-scalar mesons and the suppression 
of tensor mesons can be understood within the string model. 

Baryon production is a more complex process than the production of mesons. BB correla- 
tions in rapidity are in agreement with the assumption that the B and the B are produced 
as neighbours or next to neighbours in a string breakup. The distribution in the angle be- 
tween the baryon and the thrust axis (i.e. the general string direction) in the BB cms is 
not spherically symmetric |4], ||. Furthermore, recent data from polarized e + e~ annihila- 
tions show that baryons are more frequent in quark jets than in antiquark jets ||. Thus 
baryon-antibaryon pairs do not originate from isotropically decaying clusters. Instead the 
distribution can be understood if the B and the B are pulled in opposite directions at the 
string breakup. The simplest possibility would be that the string breaks by the production 
of a diquark-antidiquark pair, which become constituents in the baryon and the antibaryon, 
and a model based on this assumption was presented in ref J7J. In particular data on trans- 
verse momentum correlations [[| do not agree with expectations from this simple model, but 
indicate that occasionally one ore a few mesons may be produced in between the baryon and 
the antibaryon along the string. This idea was developed in the so called popcorn model ||. 
Because of the limited experimental information on baryon production, in particular con- 
cerning multi-strange baryons, some simplifying approximations were used when the model 
was implemented in the JETSET Monte Carlo [§, |9|. Thus it was e.g. assumed in the MC 
that at most one meson could be produced between the B and the B. Until recently this 
version of the MC was in fair agreement with most experimental data on baryon production. 

With the advent of new data with very high statistics from LEP it is however now possible 
to achieve a better understanding of the baryon production mechanism. In this paper we 
will present an improved model for baryon production. This model is based on the same 
general principles as the popcorn model in ref ||, but takes better account of the detailed 
kinematics of the process. The result agrees well with experimental data both for production 
rates for different baryons and for energy spectra. We also show here some predictions for a 
set of different correlations, which can be used to test the model. 

We want to mention that also other approaches to jet fragmentation have been presented, in 
particular where the production rates are mainly determined by the different hadron masses. 



This is the case for the UCLA model |T0[ , which is based on the Lund string fragmentation 
scheme, but in which the production rates are determined by the hadron masses rather than 
by the quark properties. This model gives a fair description of meson production with few 
parameters. However, it is not trivial to get a unique prediction for baryon production 
within this approach, and the most recent publication concentrates on meson production. 
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Chliapnikov et al. [I I have observed regularities between production rates and masses, 
spin and isospin for different hadrons. However, in this analysis the production of e.g. 
protons and A include decay products from heavier resonances, and it is not clear to us 
whether it is possible to find a similar relation for the directly produced hadrons. Also the 
cluster fragmentation model as implemented in the Herwig MC |12|], in which hadron rates 
are mainly determined by their masses, has difficulties to reproduce the experimental data. 
This further supports our assumption that dynamical properties deeper than just the hadron 
masses are important in the production of baryons. 

In section 2 we review the ideas of the Lund string fragmentation model and of the popcorn 
model, and also discuss its shortcomings within the default version. In section 3 we describe 
our revised model, and in section 4 we present some comparisons with experimental data. 



2 String Fragmentation and the Popcorn Model 

2.1 String fragmentation 
Longitudinal momentum distribution 

In the string model the confining field is assumed to behave like a relativistic string, i.e. like 
a vortex line in a superconductor. The model contains first a description of the decay of a 
straight string, and secondly the assumption that gluons behave as transverse excitations or 
kinks on the string. As gluon kinks have no inherent influence on particle composition, we 
will in the following discuss the fragmentation of a simple straight string. The string can 
break by the production of qq pairs, which are pulled apart by the string tension. When a 
quark meets an antiquark from a neighbouring pair, they can form a final state meson, as 
shown in Fig [I]. The production points lie around a hyperbola in x-t space, which means that 
on the average the mesons are evenly distributed in rapidity. When a qq pair is produced, 
the system is split in two causally disconnected pieces, and assuming that these pieces decay 



independently, one finds a unique solution flTj|. The probability to find a definite final state 



with meson momenta Pi is given by the following relation, if we assume that there is only 
one meson species with mass m, 

Prob. oc Y[(Nd 2 Pl 5( Pl 2 - m 2 ))5C£vk - Vtot) exp(-L4K 2 ). (1) 

i k 

Here A is the (colour coherence) area indicated in Fig [l], k is the string tension and A" and 
b are two free parameters. This expression has the form of a phase space factor (determined 
by the parameter N) times an exponential term which has the form or a Wilson area law 
(with the scale parameter b). The distribution in Eq (|I]) can be generated iteratively from 
one end of the meson system in Fig [1]. Thus, if a set of mesons are generated and each one 
takes a fraction z of the remaining light cone momentum E + p (or E — p if the system is 
generated from the other end) then z is given by the following distribution 

f( z ) = ]V (1 ~^ Q exp(-bm 2 /z). (2) 
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The parameters N, b and a are related by normalization, leaving two free parameters. With 
different quark and hadron species a larger freedom is consistent with the assumptions given 
above. Different quark flavours can be produced at different average proper times. For a 
meson with transverse mass m±, built up by flavours a and (3, the distribution in Eq @ is 
then replaced by the expression 

f a p{z) oc exp(-bm 2 Jz). (3) 

Here a is the flavour of the previously produced qq pair, and /3 the flavour of the new qq 
pair. The proper time, r, of the break-up points, or the variable T defined by 

T EE K\\ (4) 

is given by the probability distribution 

H(T) oc T aa e- bT . (5) 
From this relation we get the mean value 

(T«> = (6) 

Thus with different parameters a a , the proper time for the breakup points may depend on 
the flavour a of the qq pair produced. In fits to experimental data different a-values have 
usually only been introduced for effective diquarks in connection with baryon production. 

Flavour and transverse momentum 

Quarks with mass \i and transverse momentum k± must classically be produced in the string 
at a certain distance, d, so that the string energy between them can be transformed into the 
transverse mass of the pair 2/ij_ = 1\j[i 2 + k\ = nd. Quantum mechanically this production 
can be described as a tunneling process || |1|. The wave functions of the quark and the 
antiquark are exponentially damped in the classically forbidden region, and the production 
probability is proportional to the square of the product of these wave functions. Using the 
WKB approximation we find that this is approximately given by 

iV^ql 2 ~ exp(--/4) = exp(--/x 2 ) exp(--k 2 ± ) (7) 

Kj Kj Kj 

if the wave functions ip are normalized to f at the classical boundary. This expression 
factorizes in a flavour factor (s-quarks will be suppressed by about 1/3, while charm and 
heavier quarks can be neglected) and a Gaussian k±_ distribution. 

Vector meson suppression 

The result in Eq (0) was obtained if the wave functions were normalized to 1 at the classical 
boundary. The produced quark must however fit into the wave function of a final state 
meson. A spin-spin interaction implies e.g. that the wave function for a p-meson is more 
spread out then that of a pion, see Fig ^. Thus the wave function at the boundary is smaller 
for p, and it has been estimated that |] 

ip 2 (class, boundary) ~ l/m± (8) 
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where m± is the transverse mass of the meson. 



In this way we can understand why the p/ir ratio is suppressed from 3:1, as obtained from 
spin counting, to something closer to 1:1. The p± dependence in Eq (§) is weak for most 
mesons except for the pions, which have very small mass. For pions we expect from Eq (|8|) 
an enhancement for small p± values and a more narrow p± distribution. Due to the small 
pion mass we also expect further correlations in p± between neighbouring pions, see further 
ref 0. 



2.2 Baryon production 

As mentioned in the introduction the simplest model for baryon production would follow from 
the assumption that the string can break by the production of a diquark-antidiquark pair in a 
3-3 colour state. An antidiquark with specified flavour and spin would then appear as a new 
quark flavour in Eqs ([|[|,[7|), and different production rates would be determined by effective 
diquark masses and spin suppression factors. (Because baryons are symmetric three-quark 
states, extra weights have also to be applied in order to guarantee SU(6) symmetry.) As the 
effective masses are difficult to determine theoretically, they were represented by a set of free 
parameters in the MC implementation. These parameters determine the relative production 
ratios qq/q, udi/ud , us/ud. The MC also includes a possibility to further suppress the 
decuplet baryon states relative to the octet states, but the corresponding parameter is usually 
set equal to one. 

To improve agreement with experiments a more general framework for baryon production 
was presented in the "popcorn" model ||. In this model a diquark is produced in a stepwise 
manner and not as a single unit. Assume that a colour field is stretched between e.g. a 
red quark and an antired antiquark. The string can break if a rr qq pair is produced and 
pulled apart by the string tension. We can also imagine that a bb pair is produced as a 
virtual fluctuation. If the rb (rb) is in a colour antitriplet g (triplet g) state, the colour field 
between the produced quark and antiquark will correspond to a triplet colour field with the 
same strength as the original field. Thus equal forces in opposite directions act on the new 
quark and the new antiquark, leaving no net force on either of them. In accordance with the 
uncertainty principle they can move around freely for a time inversely proportional to their 
energy. These quarks are in the present paper called "curtain quarks" , from the picture of 
rings on a curtain-rod, sliding back and forth with no frictional losses and without changing 
the properties of the rod. In distinction to curtain quarks, the quarks that cause the string 
to break are in this paper called "vertex quarks". If the string breaks within the colour 
fluctuation region, an effective diquark-antidiquark pair is produced, see Fig [3|. 

As no net force is acting on the curtain quarks they appear as free particles. Therefore 
it was in ref || expected that the probability to find such a qq pair with transverse mass 
fi± at a separation d in space is proportional to \Ap{d, fi±)\ 2 . To produce an intermediate 
meson system with invariant mass M between the baryon and antibaryon, the distance 
d must be larger than M/k. As Aj? falls off exponentially for large d- values, a strong 
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suppression is obtained for large masses M. For this reason the possibility to have more 
than one intermediate meson was neglected in the MC implementation. The probability for 
one intermediate meson was in the program determined by the so called popcorn parameter. 
In the popcorn model the production ratios between different baryon species become slightly 
modified since a specific diquark is not always accompanied by a corresponding antidiquark. 
Thus Q production is somewhat less suppressed, as the antiquark partners of the ss diquark 
not necessarily go together in a ss antidiquark. 

Turning to the momentum distribution we note that baryon spectra in the default JETSET 
program is harder than the experimental results [15-19]. In the popcorn model the two 
quarks in an effective diquark are produced at different proper times (cf Fig|3|). The relevant 
time in Eqs (f|) and (|6|) is the latest one, the one corresponding to the "vertex quark", and it 
is natural that on average this is larger than the time for the qq production associated with 
meson production outside the colour fluctuation region. Thus we could expect a larger value 
of a for the effective diquark production (cf Eq @), which according to Eq (Q) results in a 
softer baryon spectrum. This effect is however not large enough, and in attempts to describe 
experimental data it has been suggested that first rank baryons might be suppressed [pOfl . 



(Rank denotes the order in which the hadrons are produced, counting from the string end. 
Thus in Fig [I], the meson built up by the original quark q and the antiquark q x has rank 1, 
the meson build up by qi and q 2 has rank 2 etc.) In the string model heavy baryons like 
A c or A b are only created as first rank particles (or decay products of these), and in order 
to get enough of these heavy baryons, the suppression of the first rank baryons should not 
apply to c or b jets. If the curtain quarks move along the light-cones as indicated in Fig ^|, 
a special suppression of first rank baryons may not appear natural, but, as we will show in 
the following section, it can indeed be dynamically understood in a more detailed analysis 
of the kinematics. 



3 The Revised Model 



3.1 Baryon flavours 



In the revised model we take the separate production of the curtain quarks and the vertex 
quarks more seriously. Thus we assume that the production of a vertex qq pair, which breaks 
the string, is always determined by the expression in Eqs (|?p, irrespective of whether it is 
ending up in an effective diquark or not. The s-quark suppression of all vertex quarks is 
therefore essentially determined by the K/tc ratio. 

Turning to curtain quarks we assume, following ref ||, that the production probability for 
these quarks is determined by the function A F (d, /ij_). Here is the transverse mass of the 
curtain quark and d is the minimal separation of the pair. If a meson system with invariant 
(transverse) mass M is produced between the baryon and the antibaryon, d is given by the 
sum of M and the effective diquark masses (called m\ and m.2) entering the baryon and the 
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antibaryon. Thus 

nd = M + m 1 + m 2 . (9) 

For large values of d, Ap is exponentially suppressed and as in the previous implementation 
we use the approximation 

\A F (d,n±)\ 2 ~ exp(-2/x ± d). (10) 

The relative production rates for different diquarks are thus governed by their mass dif- 
ferences. Even though diquark masses must be considered unknown, we assume that the 
differences in mass can be determined from the constituent masses, as derived in the parton 
model. 

As the quantity fj,± in Eq (|10|) is the transverse mass of the curtain quark, we note that large 
transverse momenta are strongly suppressed. Thus k± for the curtain quarks are generally 
small compared to the k± of the vertex quarks, and has therefore been neglected in the 
MC implementation. Consequently the transverse momentum of the hadron will essentially 
be given by two vertex quarks, and - like the earlier versions - the new model predicts a 
similar pj_ distribution for primary baryons and mesons []. This approximation should be 
very good for inclusive distributions, but may be somewhat too crude for correlations in p± 
or azimuthal angle. In the implementation, we let two parameters 

Az = 2 (// ±q ) /k, or (3 U and Ap = (3 8 - (3 U , (11) 

represent (fi±) in Eq fllCf ) for light and strange curtain quarks. 

From Eqs fll0|) and (]TT|) , the relative production rate for different diquarks with the same 
type of curtain quark (light or strange) can be expressed by diquark mass differences and two 
(3 parameters. The ratio between diquarks with different types of curtain quarks depends 
not only on mass differences but also on m udo . Thus e.g. 

^7 U ^°' S \ = exp(-2(/5 s m US0 - p n m udo )) 
P(ud , u) 

= exp(-2/3 s (m US0 -m udo ))exp(-2A/3m udo ). (12) 

At first this seems to call for a new free parameter. It turns out, however, that by changing 
the p s value and the decuplet suppression discussed below, one can compensate all effects of 
a variation in m u d , apart from a small variation on the A production rate. In view of the 
experimental uncertainty regarding A production, we find little need to introduce T/2 ud Q as a 
free parameter, and in the MC its value is therefore fixed to 0.2GeV. 

In the same way as for mesons, we also for baryon production expect a spin dependent factor 
related to the normalization of the hadronic wave function. Also the earlier MC contains 
one parameter for suppression of the decuplet states compared to the octet states, but in 
the default version this parameter is set to 1, because the decuplet is sufficiently suppressed 



3 Neglecting the rest mass of the quarks and using a typical value for nd w lGeV, we hnd P{k±) oc 
exp(— 2dk±) =*> (k±) ~ O.lGeV. As the fcj^-contributions to a hadron are essentially added in quadrature, 
this is very small compared to the transverse momenta of the vertex quarks, for which we have (/cj_)^ vcrtcx ^ ~ 
0.4GeV. 
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by smaller weights for production of the heavier spin 1 diquarks. In our model the three 
quarks in a baryon are produced essentially independently. Thus we expect a suppression of 
higher spin states, in a similar way for baryons as for mesons, due to the normalization of 
the total three-particle wavefunction. Following the arguments in ref QTJ] we assume that this 
factor is approximately linear in the baryon mass (cf also with the tt -7]-i]' ratios discussed 



in ref [fL4| ). Thus in addition to a suppression factor for the decuplet states compared to 
the octet states, we also expect a suppression of S compared to A due to the larger spin- 
spin interaction in ud quark pairs compared to us or ds pairs. In the MC program we have 
introduced one parameter for the decuplet /octet suppression while the extra S/A suppression 
from the mass differences (M s — Ma ~ 70GeV, M decuplet — M octet ~ 200GeV) is fixed by 
the decuplet/octet suppression through the relation 

P(decuplet) 1 P(£) 1 



P(octet) 1 + Q P(A) l + Q/3' 



(13) 



This suppression is implemented in such a way that the total production of A and £ is kept 
constant. 

Popcorn mechanism 



Using the two (3 parameters, Eq (|9|) and fllO"D, the popcorn mechanism has been implemented 



without the introduction of any new parameter. 

When the popcorn model was implemented in the JETSET MC it was assumed that the 
production of two or more mesons between the baryon and the antibaryon could be neglected, 
and there was no extra vector meson suppression in spite of the larger mass compared to 
the pion mass. This was a deliberate approximation assuming that the overestimated p and 
uj production could simulate the neglected production of two or three uncorrelated pions. 
In the new MC this approximation is not used, and the probability for a "popcorn meson 
system" is determined by the invariant mass of this system. In table [I] we show the rate of 
popcorn systems with different number of popcorn mesons for different values of (3 n . 

Popcorn systems with more than one meson could also be formed by two separate colour 
fluctuations as illustrated in Fig f|. In such a situation the baryon and antibaryon need not 
have any flavour in common. We expect this production mechanism to be suppressed, and 
it is neglected in the MC. If the production rate for each curtain quark pair is determined 
by the ratio P(qq)/P(q) ~ 0.1, the "double curtains" should be suppressed cf to the "single 
curtains" with multiple popcorn mesons by this factor. As the latter is only a small fraction 
of the total baryon production (15% for /3 U , cf tab [I]), we believe the double curtains to give a 
very small contribution, even if they are somewhat less suppressed by the factors on Eqs @ 
and (pup. In principle this kind of baryon production can be experimentally indicated if one 
can observe e.g. correlated fi~p or A ++ H + pairs, which can only be produced by the double 
curtain mechanism. 

Summarizing, apart from the s/u parameter determined by the meson rates and the diquark 
probability parameter (qq/q), the baryon multiplicities are in our model essentially deter- 
mined by three tunable parameters, @ u , A/3 and the decuplet/octet suppression, giving a 
total of five parameters. This should be compared to the earlier model and its implemen- 
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tation in JETSET, in which all together seven parameters govern baryon production (if the 
decuplet suppression is not included): The s/u and qq/q parameters plus five more, one 
suppressing strangeness in effective diquarks, one suppressing spin 1 diquarks, and three 
related to the popcorn mechanism. 



3.2 Momentum distributions 

We have assumed that the curtain quarks appear as massive particles. This implies that 
they cannot move along light-cones as in the illustration in Fig || Since the curtain quarks 
behave as free particles we assume that the curtain quark and its antiquark partner move 
along straight lines from a common starting point. We would then expect a suppression if 
some of the potential starting points lie outside the color field region. To estimate this effect 
we will be guided by a simple semi-classical model. (We use the semi-classical model mainly 
to estimate the relative suppression when the extension of the colour field is reduced. The 
overall normalization is determined by the quantum mechanical expression in Eq (fLOp.) 
As we will see, the net result is first a larger effective a-parameter for baryon production (cf 
Eq (Q)), second a suppression of first rank baryons. This suppression is stronger for light 
baryons, but rather weak for baryons containing heavy c- or b-quarks. Thus this model 



provides a dynamical interpretation of the phenomenological observations in ref [2(| . 



Let us first study a situation with only one vertex quark, i.e. a situation with no meson 
between the baryon and the antibaryon. Assume that a curtain quark- antiquark pair is 
produced in a point (called A in Fig |^a) at an invariant distance r from the vertex point C 
for the effective diquark-antidiquark pair. Fig |5|b illustrates the situation in a frame where 
the space distance between A and C is zero. (Consequently the time distance equals r.) We 
assume that the curtain quark and antiquark combine with their vertex quark partners in 
the points D and E to form an effective diquark and antidiquark. This diquark (antidiquark) 
is not massless, and thus it moves along a hyperbola with centre in the vertex point C. If 
the mass of the diquark-antidiquark pair is 2m± the distance d between D and E is given 
by d = 2m±/ k (cf Eq (|S|)). In the model we assume that the probability for this process is 
proportional to exp(— S), where S is the action. Thus, if fi± is the transverse mass of the 
curtain quarks we get a relative weight 

w = exp(— S) 

S = E*t (14) 



E = 2fi ± /Vl -v 2 = - d 2 /4t 2 . 

We note that the curtain quark and the vertex quark could also combine at a different point 
along the hyperbola, and similarly for the antiquarks. Also the two recombinations do not 
have to be at the same time in the particular frame chosen in Fig [5]b. Other possibilities 
would generally give smaller weights, and therefore we use for the general situation in Fig ||]a 
the weight in Eq fll~4] ) with t replaced by the proper time r between the points A and C. 

In a situation with more than one vertex qq pair, we have an intermediate meson system 
with (transverse) mass M. The minimum distance between the "combination points" is now 
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given by d — (M + m^i + mj_ 2 )/K (cf. Eq (^|)), and we will assume that also here the weight 
for the curtain quark production is given by Eq (|T4l) with this value of d. 

In a situation as illustrated in Fig |5|a, the total weight W is obtained by summing over 
all production points A inside the kinematically allowed region Q, defined by the condition 

r > m±/K. 

W = [ wdxdt. (15) 
Jn 

The result is a total weight W(T, p±,d) which is a function of the three variables T (the 
square of the proper time for the effective diquark vertex), p± (the transverse mass of the 
curtain quark) and d (the minimum separation for the curtain quark-antiquark pair). 

The effect of this ansatz can be anticipated: The action in Eq (|14]) becomes large both for 
large t and small t (which implies large velocities) and has a minimum for t 2 = d 2 /2. Thus 
the most favoured region for the production point A lies at some finite distance from C. The 
result in Eq ( |T0| ) should correspond to a situation when Q is large enough to contain most of 
this region, while the production probability should be further suppressed for small regions 

a 

A numerical calculation shows that for large vales of T the weight W scales with V in the 
following way 

W ~ T 5a{ ^' d) ; T large. (16) 

The power 5a(p±,d) has the effect of shifting the a-parameter in Eq (|3]) and (|5]). The 5a 
dependence on d and p± implies a unique shift of a for each diquark system. However, the 
differences are not large, and for relevant values of p± and d we find 5a m 0.5. 

The scaling result in Eq (|TjJ) is valid only when T is large compared to {hid) 2 . The weight 
W goes to zero when the area of possible production points Q is small. As seen from the 
numerical results shown in Fig || the extra suppression for small T differs for different diquark 
masses and curtain quark masses. From these differences, and the total normalization of W 
for different d and p±, production probabilities for diquarks could in principle be estimated 
within the semi-classical framework. 

In the semiclassical model the production probability goes to zero when 

n-r-w^i + in^) (17) 

goes to zero, i.e. when T goes to (nd/2) 2 . In a proper quantum mechanical treatment the 
probability should go to zero when r — » 0. T values below (nd/2) 2 should be suppressed but 
not totally excluded. Thus we feel that using the exact result in Eq ([15]) would be to stretch 
the semiclassical model too far, and instead we use for all values of and d the following 
suppression factor for small T-values 

P(T) = 1 - exp(-pr) with p w ln(2)GeV~ 2 . (18) 

This is further motivated by the phenomenological observation that the result is not very 
sensitive to the exact value of the parameter p. 
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In the case of popcorn mesons, it is not obvious how to choose the corresponding effective 
T- value to be used in the suppression factor Eq (|i8|) . In our simulations we have used the 
area indicated in Fig |] which is the largest possible area when the curtain quarks do not 
exceed the speed of light. This assumption implies that these T- values are in general much 
smaller than corresponding values for diquark vertices without popcorn, which consequently 
influences the popcorn probability. 

An important effect of the suppression for small T is a suppression of rank 1 baryons (the 
effect being stronger for light baryons than for baryons with heavy c or b quarks). For low 
rank hadrons, (T) is smaller than the result in Eq (||), which is a limiting value in the centre. 
The average values (r) for rank 1 and rank 3 baryons are shown in table |3], and for the light 
baryons (represented by p and A) we see a significant difference. This suppression of low 
rank (light) baryons obviously results in softer momentum spectra, and the results will be 
further discussed in the following section. 

For heavy baryons with a c or a b quark, the large quark mass automatically implies larger 
values for T. This is the case also if we take into account the fact that the heavy quark 
moves along a hyperbola and not along the lightcone, as illustrated in Fig |7|. (Remember 
that in the string model these heavy hadrons can only be produced as first rank particles 
including an initial heavy quark.) Thus the whole area T as defined in momentum space 
is not available for quark pair production but only the smaller area T e g, marked in the 
figure. In table |^ we have for A c and Ab presented corresponding average values, (r e g), 
and we note that these have approximately the same size as the (r)-values in the central 
region for light baryons. To obtain the values in table [3] we have been using a heavy quark 
fragmentation function which reproduces the experimental momentum spectra. This is the 



case for the fragmentation function suggested by Morris and Bowler |[21| , which is obtained 
if the space-time area indicated in Fig [7] is interpreted as a "colour coherence area" in the 
Lund string model. This fragmentation function is in very good agreement with charm and 



bottom momentum spectra [22|, and is the default option in JETSET version 7.4. Similar 



results are obtained if we use the Peterson fragmentation function fitted to the data. 

At this stage, no new tunable parameter is introduced in the model. The values of both p 
and Sa are determined by model predictions. However, the shift of the a-parameter is more 
complicated. The Sa value obtained above suppresses diquark vertices with large z values, 
i.e. vertices close to the string ends. Thus it describes the very small probability for a curtain 
quark pair to be produced around such vertices. A similar suppression is expected also for 
ordinary qq production if the pair is not massless with zero transverse momentum PJ. 
The massive qq pair is produced somewhere along two hyperbolae around the vertex and the 
range of possible production points on these hyperbolae decreases when the vertex is close 
to the string ends. This effect is taken into account when tuning the a parameter in the MC. 
Thus for baryon production there are two different effects which suppress small T-values, 
resulting in larger effective a- values. It is not at all obvious that the combined effect can 
be well described by simply adding the corresponding 5a-shifts. For this reason we have in 
our MC (as in the original JETSET MC) left Sa as a tunable parameter, with expected best 
values in the range 0-0.5. 
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4 Results 



We have implemented our model in a revised version of the JETSET Monte Carlo, and we 
will here discuss some results for e + e~ annihilation at the Z° energy. The results are also 
compared with default JETSET 7.4. 

Flavour composition 

We have made a preliminary tuning of the flavour production parameters and in table ^ we 
present the particle multiplicities obtained for /3 U = 0.6GeV -1 , A/5 = 1.2GeV _1 , decuplet/octet= 
0.19 and 5a = 0.5. We note a very good agreement between our model and experimental 
data. The discrepancy is less than 1.5 standard deviations in all cases except for E 1 * 1 (— 1.7a). 
Thus compared with JETSET default we get with fewer parameters an equally good or better 
agreement with data. 

At high energies the effective suppression of first rank baryons does not affect noticeably the 
total multiplicity of light baryons. It does however reduce the production of heavy baryons 
containing a c- or a b-quark, which in the string model can only be produced as first rank 
particles. Compared to the default JETSET version these particles are reduced to about 
70% (For A c and Ab the reduction is somewhat larger, around 60%). Due to the large 
experimental uncertainties both the larger and the smaller results are consistent with data. 

Momentum distribution 

Baryons with large momenta are suppressed for two reasons. The most important one is 
the T-suppression in Eq (|I~8|). High momenta are also suppressed by a large value for the 
parameter 5a. This suppression is, however, less effective. In Fig ^| we show momentum 
distributions for protons and A-particles for two values of 5a, and 0.5. We note that a very 
good agreement with proton and A-data is obtained for both 5a values. In Fig ^ we also 
show that the effect of T suppression is rather small on A c and Ab momentum distributions. 

Even if inclusive baryon spectra are fairly insensitive to the value of 5a, this parameter can be 
restricted if we separate the quark and antiquark jets. In Fig [TO we show baryon-antibaryon 



asymmetries compared with data from SLD ||. We note that data seem to favour a 5a 
value larger than 0, and that future studies with better statistics could be used to restrict 
the range of possible values for the 5a parameter. 

Correlations 

Correlation data can be used as a more thorough test of the model. As examples, multiplicity 
correlations for AA and AS pairs are also presented in table H The number of correlated 
AA pairs (0.106 per event) is in fair agreement with data, and the number of AS pairs 
predicted to be 0.021 per event. 



Some results for rapidity correlations are presented in Fig [11], which shows distributions in 
N = A" B g/ — |(A^bb' + ^bb') vs - \Ay\. We have compared our results with default JETSET 
and not directly with experimental data, because it is difficult to implement exactly the same 
experimental cuts. Experimental results, which have been compared with the JETSET MC, 



appear to be conflicting, as OPAL results favour a large popcorn parameter around 95% [25[ 
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while DELPHI data agree with 50% popcorn [17]. As seen in Fig 11 there are only small 
differences between the theoretical curves, but we note that to some extent the result does 
depend on the value of 5a. This is of interest because in published comparisons between 
experimental data and the JETSET MC only the popcorn parameter has been tuned. 

As mentioned above correlations in p± or azimuthal angle can be affected by the neglected 
small transverse momentum of the curtain quarks. For this reason we do not present predic- 
tions for these correlations. If wanted the p± of the curtain quarks could be implemented with 
one more parameter, which does not affect any other observable. One should note, though, 
that at higher energies, e.g. at LEP, correlations in p± or azimuthal angle are mainly deter- 
mined by jet production, while transverse momentum contributions from the fragmentation 
contributes very little. 

Energy dependence 

Due to the suppression of first rank baryons one could expect that the energy dependence 
of the baryon production might be slightly changed. This could follow from the fact that at 
lower energies a larger fraction of the hadrons are first rank. The energy dependence of the 
baryon rate depends also on the amount of popcorn, and it turns out that for the parameter 
values presented above the energy variation between lOGeV and 91 GeV is not noticeably 
different from the original default JETSET MC. 

We end by noting that our model contains two features which can be treated separately in 
the MC, first the new scheme for flavour compositions and second the suppression in Eq (|T8| ) 
leading mainly to a softer momentum spectrum. Thus in order to understand which features 
are essential for the result, it is e.g. possible to make simulations with the old model for 
flavour generation but include the suppression for small T. For heavy baryons this gives a 
smaller suppression than our full model, as the rate is reduced to ~85% of the standard 
JETSET value. 



5 Summary 

A revised version of the "pop-corn" model for baryon production is presented, in which the 
stepwise production of the quarks in a baryon is carried through consistently. As in the 
popcorn model it is assumed that a qq pair with "wrong colour" (i.e. not with the colour 
which makes the string break) can be produced as a quantum fluctuation. One or more 
secondary qq pairs can be produced within the fluctuation causing the string to break and 
producing a BB pair plus possibly also mesons. We assume here that these secondary qq 
pairs are produced in just the same way as ordinary meson-producing qq breakups in the 
string. In this way the number of free parameters can be reduced, and yet the experimental 
production ratios for different baryon species can be reproduced. 

Describing the initial "wrong colour" qq pair with a semiclassical model, we find a suppression 
of hard baryons, in particular of hard first rank baryons. This implies that the model 
reproduces well experimental momentum distributions for p and A. 
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Finally a set of predictions for BB correlations are presented. 
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relative rate of popcorn systems 


mean number of 


A, 




with N 


primary mesons 




primary mesons in 


(GeV- 1 ) 


N=0 


N=l 


N=2 


N=3 


N>3 


a popcorn system 


1.0 


0.69 


0.23 


0.06 


0.01 


0.01 


0.52 


0.6 


0.60 


0.26 


0.10 


0.03 


0.02 


0.62 


0.2 


0.46 


0.27 


0.14 


0.07 


0.05 


1.05 



Table 1: The rate of popcorn systems with different number of popcorn mesons for the same 
tune as in table [| where /3 U is 0.6GeV _1 , and similar tunes with (3 n = lGeV -1 and 0.2GeV _1 . 





Multiplicity 








JETSET7.4 


hadron 


experiments 


default 


modified 


P 


0.98±0.10 


1.20 


0.87 


A 


0.371±0.015 


0.387 


0.363 


S° 


0.071±0.015 


0.074 


0.070 


£± 


0.177±0.023 


0.141 


0.137 




0.0257±0.0014 


0.0271 


0.0274 


A++ 


0.124±0.065 


0.189 


0.118 


E(1385) ± 


0.044±0.008 


0.074 


0.056 


H(1530)° 


0.0061±0.0011 


0.0053 


0.0062 


Q- 


0.0016±0.0005 


0.0007 


0.0010 


A c 


0.075±0.024 


0.059 


0.034 


All c baryons 




0.083 


0.056 


A b 




0.033 


0.020 


All b baryons 




0.058 


0.039 


AA 


0.089±0.007 


0.090 


0.106 


AS~ 




0.019 


0.021 



Table 2: Average multiplicities from JETSET 7.4 simulations of 1.6 * 10 6 e + e events at 



Z energies, compared to experimental data, 
production model are (3 U = 0.6GeV _1 , A/3 
6a = 0.5. Experimental data are from ref [0 
of the meson parameters have been made. 



The parameter values used in the new baryon 
= 1.2GeV _1 , decuplet suppression = 0.19 and 
A c data are from ref |]26 



No extensive tune 
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hadron 


rank 1 


rank 3 


P 


1.9 


2.9 


As 


1.9 


2.9 


A c 


2.6 


(only created as rank 1) 


A b 


2.7 


(only created as rank 1) 



Table 3: (T) for light baryons and (T e ff) f or heavy baryons. JETSET 7.4 default simulation 
of e + e~ collisions at Z° energies (using Bowler correction to f(z) for heavy flavours). 




Figure 1: String fragmentation in x-t space. 
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Figure 2: Wave functions for p and tt, showing the different normalization at the classical 
turning point, which is marked with a dash. 




Figure 3: If a vertex pair q 2 q 2 is produced inside a colour fluctuation region spanned by qiq 1; 
an effective diquark-antidiquark pair has "popped out" in a stepwise manner. The model 
allows for several breakups in the colour fluctuation region, creating one or several mesons 
in between the baryon and antibaryon, as shown to the left. 
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Figure 5: Curtain pair production at point A, and a diquark vertex at point C. Figure b is 
boosted to the frame where A and C are at a common space point. 



0.5 























Mi 


m 






0.05 


0.1 






0.05 


0.2 






0.1 


0.2 






0.05 


0.5 






0.1 


0.5 






0.2 


0.5 
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Figure 6: Extra suppression P(T) = W/T Sa as a function ofT for different values of m and 
/i±. W is normed so that P fa 1 for large T. 
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Figure 7: In the case of heavy end quarks of mass /i q7 the curtain pair production area V 
cut by the hyperbolic motion of the quarks, giving rise to a smaller T e ff ^T - fi^\n(l/z). 




Figure 8: The T value for a popcorn system. 
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Figure 9: Momentum distribution -Jr- for p, A, A c and At,. x p = Phadron/Pbeam- Results 
from Monte Carlo simulations with T suppression model are shown for two different values 
of 5a, and compared to JETSET default and experiments. OPAL data are from refs [|T5 
DELPHI data are from ref [O] . 
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Figure 10: The x p distribution of 2* (B — B)/(B + B) for baryons and antibaryons in the same 
hemisphere as the outgoing quark. The model predicts more hard baryons than antibaryons 
in quark jets, as has been observed. Quantitatively, the predictions are sensitive to the 5a 
value. The plots show p and A asymmetry for 5a = and 0.5, JETSET default and SLD 
preliminary data |J. 




Figure 11: AA and pp rapidity correlations. N = N BB + iVg B — A^bb — ^bb- ^ e curves 
corresponds, as indicated, to different values of /3 U and 5a and to the traditional JETSET 
MC with default parameters and 50% popcorn. 
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